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ABS TRACT 


Variations on the traditional cost minimization of continuous review 
formulation are investigated in an effort to improve service as measured 
in terms of time-weighted shortages per unit time. It is proposed that 
the minimization of time-weighted shortages per unit time will improve 
service in current Navy Supply Operations. Various models are presented, 
without reliance upon unknown parameters such as order cost and carrying 
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cost, with necessary conditions and solution algorithms. 
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I. INTRODUCTION 


Inventories exist to provide service to customers by satisfying their 
demands from on-hand material, It follows, then, that a reasonable 
objective of inventory management is the maximization of service provided 
which is achieved by minimizing stockouts. In particular, total time- 
weighted shortages is thought to be the desired objective. 

In pursuing this objective, the manager of a realistically large, 
multi-item inventory system has a number of constraints imposed on his 
“when to buy and how much to buy" decisions. The stock points of the 
Navy Supply System have investment and reorder workload constraints which 
are real and binding. 

The classic variable cost minimization formulation is the most used 
method for solving this inventory problem, Multi-item problems are 
usually solved by assuming that they can be dealt with as a series of 
independent single item problem, In the presence of binding constraints 
on a population of items this approach is not applicable. Additionally 
the cost minimization formulation requires the estimation of cost para- 
meters which are arbitrary or at least very difficult to estimate. 

As a consequence of this argument, a series of models are formulated 
for multi-item policies subject to investment and reorder workload constraints. 
These models do not employ the standard ordering, shortage costs. This 
approach was suggested by A, P. Tully [1 | : 

In the next section, the problem formulation and the general models 
are developed. Section III develops the single item model as preparatory 


to studying multi-item case. Section IV presents a simplified multi-item 





formulation in whicn only the items reorder points are decision variables, 
The general multi-item continuous review model is developed in Section V. 


In the final section, we summarize and state some tentative conclusions, 
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II. FORMULATION 


It is desired to formulate inventory decision rules for multi-item 
inventories subject to specific constraints, The inventory decision rules 
will be of the reorder point - reorder quantity, continuous review type. 

As suggested by the introduction, the formulation to be used involves 
the minimization of total time-weighted shortages subject to: 

(1) Total average investments costs £ investment limit; and 

(2) total number of orders £ reorder workload constraints, 

Note that such a formulation would not be based on minimizing variable 
costs, 

The specific form of the model depends upon the assumption about the item 
demand characteristics and expressions for the total average on-hand inventory 
level and total number of buys per unit time. The first assumption is the 
distribution of lead time demand is normal (4, d;*) for all items, 

The following notation is used throughout the paper. For the i-th item let; 

Ci = item unit cost in dollars; 

4i = mean demand per unit time in units; 
Ali = mean lead time demand in units; 

Oi = Standard deviation of lead time demand in units; 

dcr ;) = probability that lead time demand exceeds r; 

Yr; = reorder point; and 

Qi = reorder quantity. 

Also let; 
K) = investment limit; and 


Ky = reorder workload constraint. 
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With a continuous review inventory policy an order is placed after 
the demand of Q units. It follows then that the average number of orders 
2. 


per unit time ar . For a multi-item inventory with N items, the total 


expected number of orders placed per unit time is 


‘ , 
“su (II.1) 
Tike, 

i=1 


Total inventory investment is the priced-out value of the total 
expected on-hand inventory. As shown by Hadley and Whitin | 2 | for 


continuous review the expected on-hand quantity, E(OH) is given by 
E(OH) =r +2 - +B (Qr 
(OH) 2 = ,x), 


where B(Q, r) is the expression for the expected shortages at any point 


of time. If lead time demand is normally distributed it can be shown 


(2] that 





Ber) = . [A@) - p +a) J, (II.2) 
where pir) = 5 [8* + (r-s 18 - Sef hey, 
1 - tage 
Q(x) = ve. 2 ; and 


V 27 
' @ (x) ax. 
r 


HT 


O(r) 


The expected on-hand quantity expression can be simplified by omitting 
the B(Q,r) term, and this approximation is reasonable if the risk of stock 


out is not too large. This assumption is employed throughout the thesis. 
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With this assumption, the total inventory investment is then given by the 


expression 


N 
= Ci ( ri + = ICR (II.3) 


i=] 


The expected number of backorders at any time may be explicitly 
determined from the steady state probability distribution for negative 
net inventory levels. Hadley and Whitin [2 | used this approach and 
showed that when lead time demand is normal, the time-weighted shortages 
expression is given by Eq (II.2). 

If the risk of stockout is small, then the expression for the time- 
weighted shortages can be simplified by ignoring the B(rtQ) term, which 
yields expected time-weighted units short per unit time for the ith item 


as; 
1 
re B(ri) (II.4) 


The objective can now be stated as the minimization of the time- 


weighted shortages for the entire inventory, and the formulation is: 


N 
min Z = 1 
y OL B (ri) (II.5) 
i=l 
subject to: 
N 
y ci(rit+ Se - pi) & Ks Ga 
Tell 
rs 
> 7 £% , Gi?) 
Til 


Qiy? 0, and ry;  umrestricted, 
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Note that the investment constraint will always be active, given the 
objective function used, but that reorder constraint may or may not be 


active in a given problem, 
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III. SINGLE ITEM B r) MODEL 


The basic continuous review formulation for one item with normal 


CU, Bf) lead time demand can be written as 


min: > B(r) | (ile) 
subject to: C (rt 2 -M) 2 K) ; (III.2) 
Dd ZK, OSES), 

Q 


Q> OO, and r_ unrestricted. 


A, NECESSARY CONDITIONS 
To solve the single item continuous review model, set up the 


Lagrangean function 


L(Q,r,Mm,6) = 7 B(r) + mLc(r+ 2 -)-K, +OL4 - KJ (lins4) 


OL WL 8, aL 
QQ’? Or’ O%” n 20? 


and setting the resulting expressions equal to zero yield; 


Taking the partial derivatives 


OL —- = i ALC - 62 co 

a6 92 B (r) + ; 92 0 : (III.5) 
gu = 5 Co-wg -8> (=) | +mce=o, CT) 
QL | 

mu c(r+3-f)-K =0, and (III.7) 
ge oe A a Ky = 0 B (III.8) 
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These equations may be rewritten as 


Q* . 2 [B(r) as 62_] (1.9) 
M,C 
PEA) - -m gd EH -20, (III. 10) 
c(r+3 yw) =k, (IIT.11) 
and 4. =K, . (III.12) 


These are the necessary conditions for solution. 


B. ITERATIVE SCHEME 
If the reorder constraint is active, the reorder quantity is determined 


from the equation (I.12) 


which implies: 


tt 
[Pr 


(TET 13) 


D OID 
i 
x 


~ 
No 


Note that in this case, equation (III.11) can be solved for r yielding 


Sy ee 
C 


Hence (Q*, r*) are uniquely determined from equations (III.13) and (II.14). 
If the reorder constraint is not active, the reorder quantity is 
determined along the line which is equation (III.11). It is observed that 
K 
the boundary conditions are as follows: when r = 0, Q = 2 (= + M), and 


K 
when Q=0, Lr = +Al. Plotting these conditions, Figure 1 is obtained. 
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| (Qmax, r = 0) 
h { , 
| bs - 
Kj h 


ee +h wa & ne ee een me eas iis oak 7 (Tmax, Q = 0) 


war em ae oe <n Oa om ee oe ee 
om S On cee aon aa we Oe ep ee “A we oe 


7 Ky % ' (Q*, r*) - (T max if reorder constraint 
eee Ha Ne ¥ active) 





wee = oe ow On <i <e -eee eeeh  e 
- gp oh Fee om et am aft ae, MR oe. ow oe oF 


aan Q 


K 
0 g=% Q =2(— +f) oi =— 


NO 


Figure 1. Boundary Conditicns for Single Item Case. 
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K 
It is realized that ip is greater than 2 +L), then the 
2 
investment and reorder constraints are mutually inconsistent (infeasible). 


is 


If all solutions lie on the line C(r +3 -M ) = KS for Q =e and 
2 


r&4 a hee 


2K, , then pick up a Q, solve for from the equation (III.11), 
evaluate _- and continue to iterate on Q until Bie minimum is 
found, Here parameterizing K, will give shortages as a function of the 
investment limit. 


An alternative approach would be to use a Lagrangean function with the 


boundary conditions Ez Q poet +L ) | and the equation (III.9) 
2 C 


with @=0. From these equations, it is determined that 


Wz a, 
2Ce J ay” C (lites) 


Hence a double iteration is required as the function of #,, Q, and r ; i.e., 
iterate Q,r for a given value of the multiplier. The solution is the 
minimum time-weighted shortages for the multiplier used; minimum for the 
investment level (K,) imputed to that multiplier value, ; Then a sug- 
gested procedure would be a binary search. 

Plotting shortages vs. multipliers on the Q,r and shortages planes, 


the different level of shortages will be obtained. 


C, EXAMPLE OF THE SINGLE ITEM MODEL 
Consider an item with its distribution of lead time demand normal 


(py, 0°). tet Jo= 125, @* = 25, M= 500, C = 12, K, = 720, and Ky = 8. 
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The problem is 
minimize Z= S (r) 
9 b 


subject to: c(r+3 “hy RK us 
1 


Qy0, and fr unrestricted. 


From t 4 Ko» it is determined that QZ 230 = 62.5. Hence if the 
reorder constraint is active, Q* = 62.5, r is obtained from (III.11) as 
meena and Z° = ,606. 

If the reorder constraint is not active, then a double iteration is 
required as the functionA), , Q, and © . The solution is the minimum time- 


weighted shortages for the multiplier used. Iterating on Q, and ignoring 


the reorder constraint, produces the following results: 


Time-weighted Expected number 
shortages per of units short 
unit time per unit time 


17.8446 
14,5262 * 
1692599 


20.8288 
25.0432 
32.9660 
42,6242 
53.6828 





Table 1. Time-weighted Shortages and Expected Number of Units 
Short for Single Item Case. 


Note: Two approaches seem to follow together. 
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IV. SIMPLIFIED MULTI-ITEM B Ig MODEL 


The basic continuous review formulation for the multi-item problem 


was given as: 


N 
minimize Z = ) B (ri) (1V.1) 


Qi 
i=] 
subject to: N Qi 
L 
ci(rit "> -Mi 4K, (TV.2) 
i=l 
) » ae 
a V 
a i > > (IV. 3) 
i=l 


Qi > 0, and ry unrestricted. 


Suppose the order quantities are fixed by some other criterion. 


Specifically the assumption is made that order quantities are determined 


from the equation Q; = G E ; (IV .4) 


where G is a constant, which is assumed to be the same for all items. 


From the reorder constraint (IV.3) and the assumption (IV.4) it follows 


that 
N ir, 
iCi 
C = Ko ; (IV.5) 
j=1 
or N 
> {Rici 
G — i= 
SS 


The determination of G then fixes the order quantities from equation (IV.4) 


and eliminates one set of decision variables from the problem, i.e., 
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N 


yn (IV.6) 


Ky Ci ; 


Substituting (1V.4) into equation (IV.2), it follows that 


N 


y GLUCtEist 2 ei 2. (LV.7) 
: 1 
j= : 


Equation (IV.7) can be reduced to the form 


Z 


Cit ak - + > [pica +) cifhs = K ; (IV.8) 


1 | i=1 i=l 


Mz 


Now the multi-item B(Q, r ) Model with fixed Q's can be written as 
N 

minimize Z = > Bri) 
Qi 
i=] 


subject to: 


iin 
?) 

- 

ry 

- 
iN 
4% 


i=l 
and ri unrestricted, 
N 
where Qi = . | Bic i 
Ek Nee SO | Wr 
2K. Ci 


A. NECESSARY CONDITIONS 


To solve the simplified model, set up Lagrangean function 


N , 
LC ri) = Ci ( ri) + 2 
iM : {ct 6 aly Ci ri “| (IV.9) 


i=l 


ZIM 





Taking the partial derivatives with respect to decision variables and 


setting these expressions equal to zero yields 


1 i -Mi i -/bj 
oh. = = a | ces - at) & EH Ah) i oo ) +,Cci = 0, and 


QL. _% = 
an Ci Fi K, 0. 


i=l 


Thus the conditions for solution to the problem are 


A = = | ri-pei) § ¢ ye) - Abe ery 
Ci 
OT 
ne = cH los ois aes P (rset § AS] ,  (IV.10) 
G YCLiAi 
N 
and ya es =x) (IV.11) 


B, NUMERICAL SOLUTION SCHEME 

In general, these equations cannot be solved in closed form, A numerical 
solution procedure is suggested. Let us consider a numerical method of 
solving equations (IV.10) and (IV.11). Observing equation (IV.10), which 


is 


7 = — 1 faeooe PLY (x-fi) eh 


note that the right-hand side has a lower bound of 0 since 4,30 from 


G40, ci2®iSo and fs Eh Ay : (ri-JLiy O a o O for all ri. 
Specifically} = 0 when ‘Ti becomes infinite and oO; ete - (r,;-Abi). 


d = = 0, bi: ty =(@ violates the investment constraint. 


fags 





Hence for the initial value of AL Ht iS reasonable to start wkth 


B= 0, then this implies 


= 1 
mt sar Le M) om Gc Asy| fore lei, 


< {ci ds 


or 


min | Ga 1d Pa eo) cae pg Dh. | 


Let § = pin Bas, x —— For all “1. ; 
: — se | be + LL QD ¢ 7. 
N, = § implies that there is at least one SA at its lower bound of zero. 
Hence ff] = § is a convenient starting point. Then a suggested solution 
procedure would be to begin at A) = :. solve equation (IV.10) for the 


r's and compute the value of constraint using equation (IV.12), which is 


N 
: Cc; r= HH, (IV.12) 


$ 


A bisection search will be used again, If H > K increase AL, by ae 
if H € k' » decrease N), by = . Recompute the r,'s and the value of 

constraint using equation (IV.12). If the increase (or decrease) of 
has not caused the change of inequality sign, increase (or decrease) A) 


§ 


by the same amount Ga If the sign of inequality has changed, then 
reduce the increment co and increase (or decrease) ff , solving for 
the r;‘'s at each value of ff, and computing the value of H until the 

sign of the inequality switches again. Continue until H = K,' or until 


H is within some tolerable limit of Ky . This approach will converge to 


the optimal solution rapidly. 
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From the Kuhn-Tucker theorem [ 3 | , if we have a convex objective 
function and a convex constraint region, the necessary conditions are also 
the 


sufficient. Since the constraint under consideration is linear inr , 


region is convex. To show Z(ri) is convex, consider the equation of the 


expected time-weighted shortages, 


zie 1 [at Bcri). 

Bi 
2 . 

Now if cs Zp 0 for all ri, then Zi is convex. Taking partial 
Ori 


derivatives, 





zi . 1 fe ie ri-Hi, oid -ri-Hi 
Ori (zt [cdo bah - dct] Co, 
and | 

@2 

fis 1 fo; i -Hi 

ar ey Fe gES (v.13) 


Equation (IV.13) will always be greater than or equal to zero. Under 


these conditions Ziis convex. It follows that Z is convex since it is 


the sum of convex functions. 


C, EXAMPLE OF THE SIMPLIFIED MULTI-ITEM B(Q,r ) MODEL 


Let us consider an inventory of three items, It is assumed that the 


distribution of lead time demand is normal with mean Abi and variance ig 


for the jth item, Let the item data be as follows: 


Item 1 Item 2 Item 3 

MA 1000 1500 2000 

Ci 1 10 20 

Abi ; 100 200 300 

7 100 100 200 , and let 
K, = $8,000 and Ky = 15. 


1 
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Now the problem is 


minimize Z 


subject to: 


f= 


where 


order quantities are 


With these orde:: 


minimize Z 


subject to: 


Mle 


i=l 


Notice that the 


fal 
— 
nme os eee 


G\Cciri 


3 
” 25 fires 


3 
= > ee fF (ri) 
G i Bel 
a G : 
i=l a 
= 23.6065. 


“5 


determined from equation (IV.4), 


746.5022, 
289.1190, 
236.0548 


and 


e 


the problem now can be stated as 


quantities, 


eee od 
23.6065 


- 
Il WwW 
pont 


3 
3 
Ci ri £ 8000 + S_ cifli - > tet = 11,920.5078. 
i=1 


i=l 


solution to the above problem is that vector r_ such that 


[ os apt ae peed ~ ( ri- fbi ¢ a i] . for all i, (IV.10) 
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For each value of #, , there will be some {xr} for which equation 
(IV.10) is satisfied, However, having the convex objective function, 
there exists only one {ri} for all i such that equations (IV.10) and 


8 
2 


(IV.11) are satisfied at the same time. For the initial value of M, = 


usin 1 nin __ ; - My iL _ Hi 
ing § = 2 mt pha [#te ap tbe] 


3 


for all i when ri=0, Mo = 0.0318 will be used. Decreasing M from 


0.0318, Ne = 0.0055 is obtained, which implies 


Pl =) 234.3 7505 
TQ = 264.0625, and 
¥3 = 453.1250): 


Checking the constraint, it is found that 


3 


> Ci ri = 11937.5000, 
4=1 
which is within 0.2 percent level of K) : 
The expression for the time-weighted shortage per item per unit time 


is 
Zz “+ [8 Oy + (ri fy “tg BS -1 Oj crs ey EH] 
i G 9 Ci 
The time-weighted shortage per unit time is computed as 


Ta 140)_ 100 140 | 
ee. = 1007 + (240-100 240-100 
“1 ~ 5376065 Alo00 | 3 { ) >} ace) rf 1¢G G0) 


0.2179 , 

= 10 1 2 4 ae | 
Z — — = 
2 ~ 23.6065 | 1500 [ 2 foo Beet) $83) ee 200) (oo) 


2.7634 , and 
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ae 20 [1 2 2} 15) _ 200 15 
Ee 1 + (285-300 -15 ~ 200 (285-300) (-_15 
73 33.6065 55 | 2 { 200 rie r? ( 335] 


10.5232 | 


Total time-weighted shortages per unit time for the entire inventory 


are 


5 
Zz * > Zi = 13.5045. 
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V. GENERAL MULTI-ITEM B (Q,r _) MODEL 


Let us consider the completely general continuous review model for 


multi-item B(Q, r ) Model which was: 


N 
minimize a > Wi 6 ( ri) 
—-.o ' 9 
=; 1 
subject to: : Qi 
> ci (rit 9 -Miy ex 
i=] 
 o 
i & - 
i 8 


i> 0, and ty unrestricted, 
where  B( ri) = 1% + ¢ ri- My | & EAL, - S (riding GH, 


and Wi is a weighting factor. 


A. NECESSARY CONDITIONS 


Applying the Lagrange Multiplier technique, the following results will 


be obtained: 


N wi N i . NR 
Li, rath, => Qi B ceo) C4 (ri 5 A)-x] lS +} (V.1) 
i=l i=1 


=7 


QL_-_Wi Bcriy +o . Oi g 
Q,2 2 


; (V,2) 
Q, 2 | 


er Bel crea ® (== zie =) -o: os f,] +4], Ci = 0; (V.3) 


N 
Bile = Qi - ae Ae 
2ii, y ci (ri + 47 Abi) K, = 0; and (V.4) 
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oY) 


L 


N 
a1. . ee ; 
“> ¢ =e 7” 


Then equations (V.1) and (V.2) can be solved for each multiplier yielding 





Q - [oid (ts - (ri-Miy@ A eat - 2 B( ri) 


2 Ai Wi . 
or 
Q* = o [wi Bc ri) + Os | 
i (V.6) 
M, Ci 
and : 5 ELM ~ con Jee a 
ee noe ae i (V.7) 
CiWiQi ’ 
N i 
a Ci( ri — -ALi) = K (V.4) 
i=1 


N 
ds 
> or = By ; Wes) 


Considering equation (V.7), the right-hand side of equation is always 
positive, which implies 4) 0, since 4],= 0 implies Qi = which 
violates equations (V.5) and (V.4). 

Observing equation (V.6), the right-hand side of this a is 
always positive. This suggests two possible cases for considerations. 


Case I, Ny 0, 0 = O (Ignoring the reorder workload constraint). 


The necessary conditions in this case are: 


2 2wi Bc ri) 
i" qa CO (V.8) 


29 





vf 


Ai py all (V.7) 
CiWiQi 
and N i 
I Ci Cr, +3 - Ai) “K, . (V.4) 
i=1 


Case Il. MN? D, 6>0 (Both constraints are active). 
The necessary conditions are the same as the previous page with Eq 


(V.6), Eq (V.7), Eq (V.4), and Eq (V.5). 


B. ITERATIVE SCHEME FOR CASE I 
Equations (V.8) and (V.7) can be thought of describing two curves in 
the Q-r plane. Equation (V.8) clearly shows that order quantity approaches 


zero with increasing reorder point to infinity. Furthermore, from equation 


(V.8), i? € 0 for all i; 


oe -[ oy Ay - Cri My EEA ] 
dry a nan eee. 


| /. A, ciwipC ri ra) 


7 
Taking the second partial, aa is greater than zero or less than zero 


2 
dr; 


depending upon the sign of numerator, which is 


2 
a7Qi _ 2A CiWi B(ri) § tt -At) 2 [oi cis diy - (ri- Ming =A) 


dy, 2 
[27,03 Wi B (F4) ] 





From equation (V.7), note that the nonnegative quantity 


[ “4 (ula Sit, ieee i A At) ] 


decreases to vero as ry increases to infinity, at which time the order 


quantity will be zero. However, taking the first and second partials 


20 








2 
from equation (V.7), dri < 0, ca O will be obtained; 


dQi dQ” 
a = . § cil, ra 0 fOredas 1, 
i ri- 
oo 
and _ Aciwi Oi Y 0 for alli. 


Le b (Fy we 


If one plots the two curves, Figure 2 will be obtained. 
By setting the value of rj to zero in Eq (V.7) and Eq (V.8), Qa and 
Qg will be obtained respectively, which are 
i@ ( man) (et Abi) O° ai 
A, ciwi ; 


Qa = 


and 





2. wi 1 ¢ 4? Hi] 


Ae 


To have a solution, it is necessary to have Q,€Qp. Equating Eq (V.7) 


and Eq (V.8), yields 





— Uh. rs. 2 
oa Bt — [bie - (ry -Mi) 6 mT 
Qi ye Mic (M, CiWi) 2 
or 2 
ri-Hi, BAL} 


MN, = 


ociWi> 6 (r;) 


Notice that Eq (V.9) is the function of the decision variable r only. A 


double iteration as a function off), Q andr is required, i.e, iterate 


Sy 





Eq (V.7) 


Figure 2, 


(Q*, r*) 


Eq (V.8) 


Qa 


Q as Function of r for Case I 
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Q, x for a given value of the multiplierf],. For the initial value of q , 


let yi 
i hej 72 
dip (-) +H Gc - 4] 
6 = min Loe e mace) 
for all i, pick = 2 as a convenient starting point. Then solve 


equation (V.9) for the ere} , compute {Q,'s} using Eq (V.8), and 


N Q. 
lculate th 1 ig traint using Eq (V.4). Let = Ee 
calculate e value of constraint using Eq ( ) e C; (r,t : JL ,)=H. 


i=l 
A binary search will be used. If H) Ry» increase 4 by + *. "lion @ SS ; 
decrease 4, by =. . Compute the value of H. If the increase (or decrease) 


of has not caused the sense of inequality sign, increase (or decrease) 


K, by the same amount Be: If the sign of the inequality has changed, then 
reduce the increment to a , and increase (or decrease) 4} , solving for the 
{r,'s} and f Q,'s fat each value of A] and computing the value of H until 


the sense of inequality changes again, Continue until H = “ or until His 
within the tolerable region of K)- This method will converge to the solution 
vector rapidly. 

From the Kuhn-Tucker theorem weal if we have a convex objective 
function and a convex constraint region, the necessary conditions are also 
sufficient, Since the first constraint under consideration is linear in 
r, the region is convex. To show Z(r;) is convex, it is required that 


the Hessian is greater than zero or equal to zero; i.e., 


O2i O24 
Or; Or, Or, 9G 
ZO. 
B2: ; 
1 O23 
OQ Or; OQ Q; 


Let (ry) = did — =r MG AEA) 
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Taking first and second partial with respect to each decision variable, 


the following will be obtained: 











02, = A (r;) Zs on B(r;) 
Or; Qs Qy Qi2 
Oz, | i § oe re 7 aa -_ Oz,  - 2B(ry) 
Or,Or, 2% Q, a 
Ot, _ h(i) Bt, _ A (3) 
Ori0G% 9.2 » 2 OQOr; 9,2 
From Equation (V.10), 
ri-My 
o 04; | (Oly _ Bz Oz 4G fi) 2B) _A(r,) AG,) 
Or,08, tj 29; Q. Q; 02.0 ry Ori0 Qy Qs Qe Q;2 qr ; 
i sel { (x » fe (V.11) 
a, Bap gy i 


Hence Equation (V.11) is greater than zero or less than zero depending 
upon the sign of bracketed quantity. Under these considerations, the 


convexity of Z; is unknown; only a local minimum can be assured. 


C, EXAMPLE OF THE GENERAL MODEL FOR CASE I 


Once again consider the inventory of three items from IV.C, it is 
assumed that the distribution of lead time demand is normal with mean 


2 th 
K, and variance Ox for the i item. 


Item l Item 2 Item 3 
i 1000 1500 2000 
cr 1 10 20 
Mi 100 200 300 
Pie 100 100 200, and let 


Ky = $8,000 and Ky = TS 
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Reviewing section (V.A), note that the solution vectors must satisfy 


Q? Z 2W (ra) 


M x ; (V.8) 
he ote) 
f, - 4g ci = (ry My) G : ay 
CyWiQs 


and 
N 
> Ci (ry +t fll) - =x, . (v4) 
i=1 
Using the seaeen scheme described in section (V.B), this example was 
solved with the results shown in Table 1 utilizing three investment 
levels. Comparison between the General Model for Case I and Simplified 
Multi-item Model of section III was made for the purpose of determining 
how good or bad the simplified model performed in terms of time-weighted 
shortages and marginal cost imputed to the first multiplier. These results 
are shown in Tables 2 and 3. It can be seen that the Simplified Model 
produced time-weighted shortages which were from 25% to nearly 100% larger 


than the General Model for Case I, 


D, ITERATIVE SCHEME FOR THE GENERAL MODEL 
Again equation (V.6) and V.7) can be thought of as describing two 
curves in the Q-r plane, Equation (V.6) shows that the order quantity 


approaches a fixed value, Q,; as the reorder point goes to infinity. 





Furthermore, from equation (V.6), < 0 for all i; 


i 


ay [e.g Ca ——) - (,- uppers | 


dr; 


2nc.W, [0x )+ 64% J 
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506.2500 













a 1 2 
VARIABLE 


60.8684 79.4553 178 .4906 
0.0044 0.0044 0.0044 
1.7529 7.8755 


Es reu.saae | 307.8125 | x, = 4,000 






69.4359 101.1958 247.2386 
0.0154 0.0154 0.0154 
0.5334 7.7740 37.9866 Z = 46.2940 


400.0000 687.5000 K, = 12,000 


Table 2. The General Model With an Inactive Reorder 
Workload Constraint (Case I) 








ITEM 
3 REMARKS 

VARIABLE 
264.0625 


s 234.3750 453.1250 


746.5022 237.1130 


236.0648 


hl, 0.0055 0.0055 0.0055 


3 0.2179 2.7634 10.5232 Z = 13.5045 
L 


ML 0.0173 


300.0000 343.7500 


i 
—_ 746.5022 289.1190 
7 _ 
0.4632 Z = 2.1483 


Table 3. Simplified Multi-item Model Solutions 


O 
alii 





296.8750 K 


4,000 
236.0648 kone 1s 


0.0173 





|: 


612.5000 K 


12,000 


236.0648 Ky = 15 





= 
a> 
i ; 4 7 
_ 
2a _ = - 
a => ' 
é ) 
_ 
: 
Ss | 
—-_ 
> 











d7Qi 


The second partial, “~~5 , is greater than zero or less than zero 
dr; 


depending on the sign of the numerator; 


2 


3 
dry [2qesms { Bien + ORs me 
Observing equation (V.7), note that the same properties as section 
(V.B) will be obtained, which are 


De 
ele 0. and 22! S 0 fer ell i. 


If one plots the two curves something like Figure 3 will be obtained, 
Observing Figure 3, the numerical search in Hadley and Whitin [1 | 
will be used for tie general model. To initiate the numerical procedure, 

it is necessary to determine reasonable multiplis: values. These two 
multipliers are used to compute r when Q,; = (ree. = Q. The solution 
could be started with the point where Q = a. on the curve defined by Eq 
(V.6). The r value so obtained is used in Eq (V.7) to compute a new Q 
value Qo; i.e., move from the point (Q;, rj ) on the curve defined by Eq 
(V.6) to a point on the curve defined by Eq (V.7) having the ordinate rj. 
The Q5 value is used in Eq (V.6) to compitte a new r; i.e., we move from 
the curve defined by Eq (V.7) to the curve defined by Eq (V.6) at constant 
Q. Hence a series of steps is obtained as shown in Figure 3. It is 
clear that the numerical method must converge to Q* and r* for every i, 
given Qa Z Qn 


Then compute the value of Eq (V.5) and Eq (V.4). 
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Q as Function of r for Case II. 


Figure 3. 
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In changing the values of two multipliers sequentially, 4}, = G20, 
tno2, .,.,andmon= 152, ...,aeset of solution vectors, {Q, rf} will 
be obtained, ‘Then compute the time-weighted shortages and the value of 
Eq (V.5) and Eq (V.4). 'Good' solution vectors can be found with small 
variance of Eq (V.5) and Eq (V.4) and approximate the minimum time- 
weighted shortages, Once obtaining this, iterate with more careful 
selection of multipliers, compute the objective function and check the 
variances of Eq (V.5) and Eq (V.4). Continue iterating on the multipliers 
until tolerable dipeeeentes basween the right and left sides of Eq (V.5) 
and Eq (V.4). 

More sophisticated approaches are available. Among these is a search 
technique proposed by Fiacco and McCormick [4] , and a two state variables 
optimal allocation using dynamic programming. 

The sequential unconstrained minimization technique [ 4 | , oUME;, is 


N 
based on the minimization of a new function P(X,?) = £(X) “th 1/g; (X) 


i=l] 


over a strictly monotonic decreasing sequence of C- values {Px} ; 
Under certain restrictions that will be reviewed subsequently, there 
exists a sequence of feasible points {xh )} that respectively minimize 
frcx ; pt , and it follows wee x ¢ Pi=> X , a solution of original 
function as Qe. —> 0 (k->0 ). The following is a concise summary of 
the steps describing the computational algorithm: 

1, Select a point X° interior to the feasible region, 


2. Select ae the initial value of P using 


Pr - Fe) gy wR) / | ype], 
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or 


de | 
SJE(K )Ha 9 £(X°) 1 


-1 
V(x) ny p(x?) 


where p(X) = > 7 ifsc, H, and Ho are the Hessians of f(X) and 


p(X) respectively. 


3. Determine the minimum of P(X, Pe) for the current value of P,. 


using Gradient Methods, 


1 
x - OVP (x), 


ms 
ii 


or 


oa g(r fx, x5] Vv P(x’) 


4. If k>{ , estimate solution using extrapolation formula. 
5. Terminate computations if final convergence criteria are satisfied, 
£(X) ~ c[x@), ho ]K< € . The theoretical optimum value Vo is bounded 


by the dual and primal function values respectively, 


of xp, Acer] Sv. eLxcey] . 


If the value v, is not within the bounds above, go to step (6). 


6. Select ec = Pfc , where C> 1. 
7. If k>} 1, estimate minimum for reduced ( - value, using an extra- 


polation formula, 


E. EXAMPLE OF THE GPNERAL MODEL (CASE IT) 


Consider an inventory of three items with the following characteristics: 


Item l Item 2 Item 3 
4 1000 1500 2000 
Ci 1 10 20 
My 100 200 | 300 
032 100 100 200 
K, = $8,000 Ky = 15 


AV 





The search technique proposed by Fiacco and McCormick | 4 | was used, 


The following initial feasible solution was used: Q = 600, r, = 200, 
= = = = — 3 
Q2 = 270, r, = 260, Q, 300, r., 400. Zi B(ri) = 17.8084 and 
i=1 Qt 


f 


- 1,6991 were obtained, 


This example was solved with the following results, 


Q, = 362.718, Dg 22-1433), 
Q, = 266.403, ry = 272.0930, 
Q, = 302.6870, r, = 425.7958, P, = - 0.00002231, 


and the value of the objective function (Z) was Z = 13.39781. 
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VI. SUMMARY AND CONCLUSIONS 


It is apparent that the general models proposed do not represent the 
ultimate answer in multi-item inventory theory. Since the conditions of 
concavity and convex sets are not maintained for the general model, the 
computational algorithms which are based on concavity and convex sets of 
objective and constraing respectively could not apply directly. Though 
the P function in section (V.D) Bonen prohibitively difficult to work 
with computationally, it has been understood that it is minimized efficiently 
and accurately for the great preponderance of problems solved to date, by 
means of the second order optimum gradient method. 

From the example problem we see the general problem in section (V.E) 
provides solution which is better than the solution of the simplified 
continuous review model in section (IV.C). The major advantage of the 
simplified continuous review model is its computational ease, 

While an efficient algorithm for solution of the general problem has 
not been presented, the advent of high speed computer has opened this field 
of numerical iterative procedures for large inventory systems. 

Since the value of the constraints are more easily determined than the 
order cost and holding cost, the model proposed seems much more appropriate 
than the traditional variable cost minimization models. 

From the example problem given in section (III), it is apparent that 
time-weighted shortages and expected number of units short per unit time 
seem to follow each other closely, 

Finally, all the solutions to example problems were obtained by 


utilizing a digital computer. 
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